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O ' Id th e no-core shell model formalism we compute effective one- and two-body operators, 

^ , using the Lee-Suzuki procedure within the two-body cluster approximation. We evaluate 

£Z< ■ the validity of the latter through calculations in reduced model spaces. In particular, 

we test the results for the two-body system and find that indeed the effective operators 
in the reduced space reproduce the expectation values or transition strengths computed 
, in the full space. On the other hand, the renormalization for operators in the case of 

% j ■ 6 Li is very weak, suggesting the need for higher-body clusters in computing the effective 

' interaction. 



1. Introduction 

One of the important steps in validating any theoretical model is a good descrip- 
tion of the experimental data. This has to include, aside from the energy spectrum, 
observables and transition strengths, which provide an important test of the theo- 
retical wave functions. 

Starting from first principles, the no-core shell model (NCSM) has been suc- 
cessful in reproducing with very good accuracy the observed energy spectrum in 
light nuclei. Thus, within the NCSM framework, one starts with a realistic nucleon- 
nucleon (NN) interaction, 1 more recently adding also three-nucleon forces, 2,3 and, 
by employing a similarity transformation, briefly reviewed in Sec. 2, obtains an 
effective interaction in a restricted model space. The initially infinite dimensional 
problem then becomes numerically tractable, so that diagonalization in a large, 
translationally invariant basis provides the theoretical spectrum as well as the wave 
functions of ground and excited states. If the transformation is exact, the energy 
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spectrum and expectation values of other observables are preserved. 

Much effort has been directed toward computing effective interactions, and with 
increasing computing power and more efficient shell-model codes, we are able nowa- 
days to describe heavier and heavier systems in larger and larger model spaces. 
Less progress has been made, however, in the direction of consistently computing 
other effective operators within the NCSM. The only operators thus computed have 
been the effective point-proton radius, 1 relative kinetic energy, 4 and the NN pair 
density 4,5 However, important operators such as semi-leptonic weak and electro- 
magnetic have not been consistently treated in the NCSM, and this is a possible 
explanation for the failure to describe adequately transition strengths. 1 

Historically, the method to obtain agreement with the experimental transition 
strengths was to use bare operators, with enhanced (e.g., E2) or quenched (e.g., 
Gamow- Teller) effective charges or couplings instead of the bare ones. But pertur- 
bation theory fails to predict the needed phenomenological charges; 6 considerable 
success was however reported within the NCSM framework in a restricted model, 
where the effective quadrupole charges have been shown to have the expected values 
in a single harmonic oscillator shell. 7 In the present paper, we test an implementa- 
tion of general effective operators within the NCSM, by performing calculations in 
two- and many-body systems in limited model spaces. 



2. Formalism 



In the NCSM, all nucleons are treated on an equal footing. In this work, we consider 
two-body interactions only, that is, the intrinsic Hamiltonian describing the many- 
body system is 



Ea = tX 

i>j 
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NN 
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i>3 



where m is the nuclcon mass, and V£j the bare NN interaction, such as the 
Argonne potentials in coordinate space 8 or the non-local CD-Bonn. 9 

Adding a center-of-mass (CM) Hamiltonian to the ^4-body Hamiltonian (1) 
does not change the intrinsic properties of the system. We bind the CM through a 
harmonic oscillator (HO), so that the total Hamiltonian can be cast in the form 
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This introduces a pseudo-dependence on the HO frequency O; however, the cluster 
approximation described below will introduce a real dependence on f2. 

Following Da Providencia and Shakin 10 and Lee, Suzuki and Okamoto, 11 wc 
construct a similarity transformation able to accommodate the short-range two- 
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body correlations by introducing an antihcrmitian operator S, so that the effective 
Hamiltonian H is given by 

H = e~ s H^e s . (3) 

In principle, both S and H are ^4-body operators (even if originally we started 
with only two-body interactions), but determining S exactly would be as hard as 
solving the original problem. We will return to the calculation of S below; for now 
we assume that no approximation has been made. 

The infinite Hilbcrt space associated with the system can be split into the finite 
model, or P-space, and the complementary, or Q-space, with the projectors P and 
Q spanning the entire space, P + Q = 1. We determine the previous transformation 
by requiring that the new Hamiltonian completely decouples the P and Q spaces, 

QHP = 0, (4) 

in addition to the conditions 12 PSP = QSQ = 0. With these restrictions, the 
operator S can be formally written 12 by means of another operator uj as 

S = arctanh(w — u;t), (5) 

with QojP = u. In terms of the new operator w, the effective Hamiltonian in the P 
space becomes 

H eff-PHP- =H A r——^ =, (6) 

and, analogously, any arbitrary operator can be written in the P space as 

V P + VP + u'uj 

In order to determine the effective operators, one has to compute the transfor- 
mation operator lu, which is equivalent to determining S. But before addressing 
this problem, we would like to comment on the last two equations. We mentioned 
in the introduction that very little progress has been made regarding effective op- 
erators, while one can compute with very good accuracy the effective interactions. 
While it is true that Eqs. (6) and (7) are formally similar, the difference is that 
the interaction is given in relative coordinates, while a general operator is usually 
defined in the single-particle basis. For the effective interaction (6), the problem can 
be broken down and solved one T) channel at a time, the dimensions involved 
being much smaller than for the initial problem. The same procedure cannot be 
applied, however, for a general operator in the single-particle basis. The solution is 
to transform ui to the single-particle basis, but this involves a complicated approach 
which includes a large number of states. The problem simplifies again for special 
operators which are given in relative coordinates, such as the point-proton radius, 
relative kinetic energy, or NN pair density, and these have been investigated, 1 ' 4 ' 5 
as previously noted. Finally, note that because of Eq. (7), a one-body operator will 
end up having non-zero genuine two-body elements in the model space. 
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We return now to the problem of determining u>. Note that this operator con- 
nects vectors in P to vectors in the Q space. Therefore, a simple way to compute 
u) is 1 

(a Q \(j\a P ) = ^2(a Q \k)(k\a P ), (8) 
k 

where \otp) and \<xq) arc the basis states of the P and Q spaces, respectively, 
the |fc) are the eigenvectors of the Hamiltonian in the full space H^\k) — E/~\k), 
and (ctp\k) is the matrix element of the inverse overlap matrix (cep\k), that is 

E Qp ( fc '|ap>(«p|fc> = <W- 

Note that computing uo by means of Eq. (8) requires knowledge of the solutions 
to the initial ^4-body problem, that is the eigevectors \k), which are actually our 
goal. While formally everything was exact up to now, in order to use this procedure, 
we have to make one approximation. That is, instead of using the eigenvectors \k) 
for the ^4-body problem, we solve these for the a-body system, with a < A. This 
is called the cluster approximation, and in the limit a — > A, the solution becomes 
exact. We point out that the cluster approximation induces a real dependence of 
the effective interaction upon the center of mass HO frequency; the solution to this 
problem is to search for a range of values over which the results vary only slightly, 
i.e., are weakly Q dependent. 

The interaction obtained by means of Eq. (6) at the a-body level is then diago- 
nalized by means of the Lanczos algorithm in the many-body space spanned by an 
anti-symmetric, translationally invariant basis of A particles. As noted, if a = A, 
one obtains the exact solutions. The way to reach convergence to the full space 
solution for a < A is by increasing the dimension of the model space P, keeping 
the cluster approximation at the same level. (Note that when P approaches the full 
space, the effective interaction approaches the bare one.) Although results calcu- 
lated at the three-body cluster level have been recently reported, 3 ' 14 in this work 
we restrict ourselves to the two-body cluster, or in other words, a — 2. (Three-body 
cluster calculations are computationally very demanding for obtaining the nuclear 
spectra and much more so in the case of other operators.) 

Finally, when computing the effective operators by means of (7), we have to make 
one additional approximation. That is because, as mentioned before, for arbitrary 
operators one cannot work in the relative system, and one needs to go to two-body 
matrix elements described in the single particle basis. This involves a large number 
of two-body matrix elements from the Q space, and even for the simplest case one 
cannot take all of them into account. The solution is to include a restricted number 
of basis states from the Q space, that is one HO shell at a time, and observe the 
convergence in the values of interest. For example, for a 2MI P-space calculation, 
we include only states from 6fi.fi, etc. In each case we compute the two-body 
matrix elements of the effective operator, and when we do not see a change in their 
values, we assume that we have reached convergence. 
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Fig. 1. Transition strengths of effective E2 and Ml operators, and the ground-state expectation 
value of the effective kinetic energy: evolution with the size of the Q-space included. The left-most 
values were obtained utilizing bare operators in the 2hQ model space, while the straight lines mark 
the exact value obtained in the full space. 



3. Results and discussion 

We apply first the procedure described above to the two-body system. In this case, 
the two-body cluster is the exact transformation, obtaining in the restricted space 
the exact expectation values or transition strengths as in the full space. 

Note that w is a two-body operator, and, therefore, in order to compute effective 
one-body operators, we have to write the latter in a particle-number dependent two- 
body form: 

0« = 5>|0«|6>ata 6 = Y.^\0 (1) \b)5 cd a\a\a d a b , (9) 

ab abed 

with the appropriate anti-symmetrization of the two-body matrix elements. In Eq. 
(9), the indexes of summation run over all single particle states, a* and a are 
the usual fermion creation and annihilation operators, respectively, and 8 is the 
Kronecker symbol. 

In the two-body system, the two-body cluster approximation is exact; conse- 
quently, we have to regain in the restricted space exactly the values in the full 
space. For testing purposes, we assume that the full space is the lOfuo shell-model 
space, in which we use the bare, isospin-independent Argonne V8' interaction. 8 We 
emphasize that although we specifically employ one proton and one neutron in our 
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Table 1. °Li: B values of select E2 and Ml transitions, and expectation values of bare and effective 
relative kinetic operators. For comparison, we list the corresponding results with bare operators 
in the full space. 



Model Space 


2ftfi 


AMI 


Full Space 
(8M7) 


Bare Op. Eff. Op. 


Bare Op. Eff. Op. 


£(£2; 1+0-^ 1+0) 
B(E2; 3+0 -» 1+0) 
B(E2- 1+0 - 1+0) 
B(E2;2+0 -» 0+1) 


0.0080 0.0079 
1.9512 1.9315 
0.8855 0.8765 
0.0870 0.0866 


0.0097 0.0097 
1.7575 1.7548 
0.7130 0.7117 
0.0865 0.0867 


0.0364 
1.5750 
0.5791 
0.0820 


B(M1;0+1 -» 1+0) 


5.3926 5.3549 


5.4209 5.4051 


5.3346 


(l+0|T re; |l+0) 
(0+l|T re( |0+l) 


116.15 118.46 
116.38 118.58 


124.66 125.84 
124.56 125.67 


135.84 
135.19 



calculations, this space does not give a realistic description of the deuteron, and is 
for testing purposes only. We then transform to a 2hw model space, computing the 
effective interaction. Following the procedure described in the previous section, we 
compute E2 and Ml transition operators, as well as the effective relative kinetic 
energy by including one shell at a time from the Q space. Figure 1 shows that 
indeed when one includes most of the states in the Q space, the effective operators 
reproduce the values in the full space. One might argue that the convergence is 
very slow, and that one needs most of the states to reach convergence. This is, 
however, most likely an artifact of the smallness of the full space, which makes the 
contribution of most Q-states equally important. 

In many-body systems, the situation is different from the two-body system be- 
cause the transformation from the full to the model space is not exact. We in- 
vestigate the influence of the two-body cluster approximation in 6 Li, by assuming 
that the full space is 8Mj and the interaction has a Gaussian form (due to the 
short-range repulsive core of realistic bare potentials, we cannot perform meaning- 
ful calculations in such a small space). We then compute the effective interaction 
and operators in 2 Ml and AMI model spaces, and summarize the results in Table 1. 
Note that because the initial full space is rather small, we have a more important 
contribution from higher order clusters. We mention however that although the 
ground state energy for the AMI model space still differs by a few percents from the 
full space (8M1), the spacing between the states is about the same. All shell-model 
calculations have been performed with a descendant of the MFD code. 15 

Table 1 shows insignificant differences between the results obtained with bare or 
effective operators (note that in this case we show only the values when we included 
the full Q space). We have to point out that previous calculations have produced 
very strong renormalization of the kinetic energy operator in 4 Hc 4 and 12 C, 16 but 
here we find it to be very weak. The explanation lies again in the smallness of the full 
space which makes the higher-order clusters, most likely up to six particles, much 
more important than in realistic calculations. Nonetheless, we have unpublished re- 
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suits in realistic calculations in 12 C for effective electromagnetic operators, 16 which 
show similar behavior as the results in Table 1. This suggests that for such operators 
the two-body cluster does not produce enough renormalization, to be confirmed by 
further investigations in realistic model spaces, using realistic interactions. Note 
that the relative kinetic energy is a short range operator, while the E2 transition 
and point-proton radius (for which was also found a very weak renormalization) 
operators are long-range operators. It is, however, known from effective interaction 
calculations that the two-body clusters renormalize the short-range correlations of 
the NN interaction, so it follows that the renormalization for long-range operators 
should be weaker. What is somehow surprising is the extremely small renormaliza- 
tion of the E2 operator, which should have also a short-range component (there is 
no surprise that the Ml operator does not renormalize much, as it does not mix 
different shells, and, hence, the contribution from outside the model space should 
be very small). 

4. Conclusions 

We have implemented the NCSM formalism to general one- and two-body opera- 
tors. For effective interactions we use a procedure which is equivalent with using 
the entire Q space (for details we refer the interested reader to previous work 1 ). 
Arbitrary operators present a significant difference, as, in general, one cannot use 
the whole complementary Q space. Hence, we developed a convergence procedure 
by adding Q states one shell at a time. While in the calculation presented the con- 
vergence appears to be slow, we expect that this is due to the limitation of the size 
of the model space, where higher-order clusters play an unrealistically big role. 

In the two-body system, where there is no approximation in computing the 
transformation from the full to the model space, we regain in the exact B(El), 
B(Ml), or expectation values by employing effective operators. Furthermore, these 
renormalized results differ significantly from the ones obtained by using only the 
bare operators. 

Our test in a many-body system, where we employ only the two-body cluster 
approximation, shows a weak renormalization for electromagnetic operators as well 
as for the relative kinetic energy. The latter has been previously computed in real- 
istic systems, and the renormalization was found to be significant. We explain the 
present results by the size of the full space, which is very small in our example, 
making the higher-order cluster contributions more important than in realistic sit- 
uations. However, for electromagnetic transition operators, we found similar results 
in a realistic calculation for 12 C, 16 which would suggest that, compared with the 
effective interaction, the higher-order clusters play a bigger role in computing long- 
range effective operators. This is a hypothesis remaining to be verified in realistic 
calculations. As a general caveat, any truncation of the space could induce effective 
operators with non-negligible higher-body correlations. Consequently, the renormal- 
ization properties of each operator, especially those with little or no experimental 
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data available (e.g., double-/? decays), need to be studied separately. 

Nevertheless, the theory of effective operators is important for at least two sets 
of next generation experiments. The first is the investigation of nuclei far-from- 
stability through the projected rare isotope accelerator (RIA), with important ap- 
plications to nucleosynthesis. The second is to double-/3 decay experiments, which 
will provide data on the mass of the lightest neutrino. 17 This motivates us even 
more to pursue this issue further, by investigating, for example, semi-leptonic op- 
erators at finite momentum transfer which should have short-range components, 
renormalizable at the two-body cluster level. 
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